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ABSTRACT 

Context. The cross correlation method (hereafter CC) is widely used to derive the radial velocity curve of Cepheids when the 
signal to noise of the spectra is low. However, if it is used with the "wrong" projection factor, it might introduce some biases in 
the Baade-Wesselink (hereafter BW) methods of determining the distance of Cepheids. In addition, it might affect the average 
value of the radial velocity curve (or 7-velocity) important for Galactic structure studies. 

Aims. We aim to derive a period-projection factor relation (hereafter Pp) appropriate to be used together with the CC method. 
Moreover, we investigate whether the CC method can explain the misunderstood previous calculation of the "K-term" of 
Cepheids. 

Methods. We observed eight galactic Cepheids with the HARPS** spectrograph. For each star, we derive an interpolated CC 
radial velocity curve using the HARPS pipeline. The amplitudes of these curves are used to determine the correction to be 
applied to the semi-theoretical projection factor derived in Nardetto et al. (2007). Their average value (or 7-velocity) are also 
compared to the center-of-mass velocities derived in Nardetto et al. (2008). 

Results. The correction in amplitudes allows us to derive a new Pp relation: p = [—0.08 ± 0.05] log P + [1.31 ± 0.06]. We also find 
a negligible wavelength dependence (over the optical range) of the Pp relation. We finally show that the 7-velocity derived from 
the CC method is systematically blue-shifted by about 1.0 ± 0.2km s _1 compared to the center-of-mass velocity of the star. An 
additional blue-shift of 1.0km s _1 is thus needed to totally explain the previous calculation of the "K-term" of Cepheids (around 
2km s" 1 ). 

Conclusions. The new Pp relation we derived is a solid tool for the distance scale calibration, and especially to derive the distance 
of LMC Cepheids with the infrared surface brightness technique. Further studies should be devoted to determine the impact of 
the signal to noise ratio, the spectral resolution, and the metallicity on the Pp relation. 

Key words. Techniques: spectroscopic - Stars: atmospheres - Stars: oscillations (including pulsations) - (Stars: variables): 
Cepheids - Stars: distances 



1. Introduction 

The Baade-Wesselink (hereafter BW) method of deter- 
mining the distance of Cepheids was recently used to cal- 

* Based on observations made with ESO telescopes at the 
Silla Paranal Observatory under programme IDs 072.D-0419 
and 073.D-0136 

** High Accuracy Radial velocity Planetary Search project 
developed by the European Southern Observatory 



ibrate the period-luminosity (PL) of Galactic Cepheids 
(Fouque et al. 2007). The basic principle of this method 
is to compare the linear and angular size variation of a 
pulsating star in order to derive its distance through a 
simple division. The angular diameter is either derived by 
interferometry (for e.g. Kervella et al. 2004, Davis et al. 
2008) or using the infrared surface brightness (hereafter 
IRSB) relation (Gieren et al. 1998, 2005a). However, when 
determining the linear radius variation of the Cepheid by 
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spectroscopy, one has to use a conversion projection factor 
from radial to pulsation velocity. This quantity has been 
studied using hydrodynamic calculations by Sabbey et al. 
(1996), and more recently Nardetto et al. (2004, 2007). 

Following the work of Burki et al. (1982), we showed 
in Nardetto et al. (2006, hereafter Paper I) that the first 
moment of the spectral line is the only method which 
is independent of the spectral line width (average value 
and variation) and the rotation velocity of the star. The 
centroid radial velocity (RV C ), or the first moment of the 
spectral line profile, is defined as 



RV C 



/line ^(X)dX 



(i) 



We thus used this definition of the radial velocity in paper 
two of this series (Nardetto et al. 2007, hereafter Paper II), 
to derive a semi-theoretical period-projection factor (here- 
after Pp) relation based on spectroscopic measurements 
with the HARPS high resolution spectrograph. This rela- 
tion was derived from the specific Fel 4896.439 A spec- 
tral line which has a relatively low depth for all stars at 
all pulsation phase (around 8% of the continuum). It was 
shown that such low-depth is suitable to reduce the uncer- 
tainty on the projection factor due to the velocity gradient 
which takes place between the photosphere (corresponding 
to angular diameter measurements) and the line-forming 
region (corresponding to the radius estimation from spec- 
troscopic measurements). 

In the cross-correlation method (hereafter CC method) 
a mask (composed of hundreds or thousands) of spectral 
lines is convolved to the observed spectrum. The result- 
ing average profile is then fitted by a Gaussian. In such a 
method, there is first a mix of different spectral lines form- 
ing at different levels (more or less sensitive to a velocity 
gradient). Second, the resulting velocity can be dependent 
of the abundances or effective temperature (through the 
line width), or the rotation of the stars. Third, in paper 
III of this series (Nardetto et al. 2008) we derived cali- 
brated center-of-mass velocities of the stars of our HARPS 
sample. By comparing these so-called 7-velocities with 
the ones found in the literature (generally based on the 
CC method) and in particular in the Galactic Cepheid 
Database (Fernie et al. 1995), we obtained an average cor- 
rection of 1.8 ± 0.2 km s _1 . This result shows that the 
"K-term" of Cepheids stems from an intrinsic property of 
Cepheids. But, it docs show also that the cross-correlation 
might introduce a bias (up to a few kilometers per second) 
on the average value of the radial velocity curve. 

After a careful definition of the projection factor (Sect. 
[2|) , we apply the cross-correlation method to the Cepheids 
of our HARPS sample (sect. [3]), in order to derive a period- 
projection factor relation appropriated to the CC method 
(Sect.|4|). As the HARPS pipeline also provides the cross- 
correlated radial velocities for each spectral order, we take 
the opportunity to study the wavelength dependence of 
the projection factor law (Sect.©- And finally, we quantify 
the impact of the CC method on the 7-velocities (Sect. [6]). 



2. Definition of the "CC projection factor" 

In this section, we recall some results obtained in paper II 
and we define the projection factor well suitable for the 
cross-correlation method. In paper II, we defined the pro- 
jection factor as: 



P = 



ARVc 



(2) 



where AV° is the amplitude of the pulsation velocity curve 
associated to the photosphere of the star. ARV C is the 
amplitude of the radial velocity curve obtained from the 
first moment of the spectral line. Because of the atmo- 
spheric velocity gradient, ARV C depends on the spectral 
line considered. Using a selection of 17 spectral lines, we 
thus derived an interpolated relation between ARV C and 
D, where D is the line depth corresponding to the mini- 
mum radius of the star : 



Ai?V c = a D 



(3) 



This relation was then used to quantify the correc- 
tion (/grad) to be applied on the projection factor due 
to the velocity gradient (see Eq. 3 of paper II). The 
Fe I 4896.439 A spectral line (which forms close to the pho- 
tosphere) was found to provide the lowest correction. The 
amplitude of the radial velocity curve corresponding to the 
Fel 4896.439 A spectral line was finally used (see / gra d m 
Tab. 5 of paper II) to derive the semi-theoretical Pp rela- 
tion. It is defined (Eq.© as AJ?F C [4896] = a D i896 + b , 
where ao and foo are indicated in Tab. 3 of Paper II. -D4896 
is derived from the interpolation of the line depth curve at 
the particular phase corresponding to the minimum radius 
of the star (i.e. when RV C corrected from the 7-velocity 
~ 0). L> 4896 and ARV C [4&96] arc given in Tab. Q] of this 
paper. 



Table 1. The cepheids studied listed with increasing pe- 
riod. The depth of the Fe I 4896.439 A spectral line and 
the corresponding value of the amplitude of the radial ve- 
locity curve (see Eq. [3]) are given. 



Cepheid 




I?4896 


ARVcime] 




[days] 


[%] 


[km s _1 ] 


R TrA 


3.38925 


6 


28.6±o.5 


S Cm 


4.68976 


6 


33.5±o.5 


YSgr 


5.77338 


5 


34.3 ±0 .5 


(5 Dor 


9.84262 


10 


31.7±o.5 


C Gem 


10.14960 


14 


25.5 ±0 .5 


RZ Vel 


20.40020 


4 


47.6 ±0 .5 


£ Car 


35.55134 


13 


32.8±o.5 


RS Pup 


41.51500 


4 


42.4 ±0 .5 



The corresponding Julian dates (T ) can be found in 
Paper II. 
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The projection factor suitable to the cross-correlation 
method (hereafter p cc ) is then simply : 

Ai?l/ C [4896] 



Pc 



p- 



Pfc 



(4) 



ARVcc 

where ARV CC is the amplitude of the radial velocity curve 
obtained with the cross-correlation method, and / cc the 
correcting factor to be applied. Our definition of p cc is 
independent of the 7-velocities. 

3. The CC method applied to HARPS 
observations 

We consider eight Cepheids which have been observed 
with the HARPS spectrometer (R = 120000): R Tra, 
S Cru, Y Sgr, /? Dor, C Gem, RZ Vel, I Car, RS Pup. 
Information about observations (number of measure- 
ments, pulsation phases) can be found in Paper I. 

We apply the HARPS pipeline to our data in order to 
calculate the cross-correlated radial velocities (Baranne et 
al. 1996, Pepe et al. 2002). The basic principle of the CC 
method is to build a mask, made of zero and non-zero 
value-zones, where the non-zero zones correspond to the 
theoretical positions and widths of thousands of metal- 
lic spectral lines at zero velocity, carefully selected from 
a synthetic spectrum of a G2 star. A relative weight is 
considered for each spectral line according to its depth 
(derived directly from observations of a G2 type star). An 
average spectral line profile is finally constructed by shift- 
ing the mask as a function of the Doppler velocity. The 
corresponding radial velocity is derived applying a classi- 
cal x 2 minimization algorithm between the observed line 
profile and a Gaussian function. The whole profil is con- 
sidered in the fitting procedure, not only the line core. 
The average value of the fitted Gaussian corresponds to 
the cross-correlated radial velocity (hereafter RV CC ). The 
HARPS instrument has 72 spectral orders. The pipeline 
provides RV CC averaged over the 72 spectral orders, or in- 
dependently for each order. We first use the averaged val- 
ues and the corresponding uncertainties. 

The RV CC curves are then carefully interpolated using 
a periodic cubic spline function. This function is calcu- 
lated either directly on the observational points or us- 
ing arbitrary pivot points. In the latter classical 
minimization process between observations and the in- 
terpolated curve is used to optimize the position of the 
pivot points (Merand et al. 2005). For Y Sgr and RS Pup, 
pivot points are used due to an inadequate phase coverage. 
When the phase coverage is good (which is the case for all 
other stars), the two methods arc equivalent (Fig.[lJ. From 
these curves we are finally able to calculate ARV CC (Tab. 
[2]) . The statistical uncertainty on ARV CC is set as the aver- 
age value of the uncertainty obtained for all measurements 
over a pulsation cycle of the star. 

4. A Pp relation dedicated to the CC method 

From Ai?14[4896] and ARV CC obtained for all stars we 
derive the correcting factor / cc using Eq. [2j The result is 
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Fig. 2. (a) The correction factor f cc induced on the pro- 
jection factor by the cross-correlation method is shown as 
a function of the logarithm of the period of the star, (b) 
The period-projection factor (p) relation from paper II 
(crosses and solid line) and the corrected relation suitable 
for the cross-correlation method p cc = pf cc (diamonds and 
dashed line). 



plotted as a function of the period in Fig. [2^. No partic- 
ular trend is found. However, the f cc correction factors 
are clearly statistically dispersed around a mean value of 
0.93 ±0.02. 

Following our definition (p cc = pf cc ), the corrected 
projection factors suitable for the CC method are given in 
Tab. [2] The relation between the period and p cc remains 
clear according to the statistical uncertainties: 



p cc = pfcc = [-0.08 ± 0.05] logP + [1.31 ± 0.06]. (5) 

The corresponding reduced x 2 is 1.2. We refer to this 
relation in the following using Pp cc . We recall that the Pp 
relation we found in Paper II dedicated to the Fel 4896 
spectral line was : p = [-0.064 ± 0.020] logP + [1.376 ± 
0.023]. These two relations are shown in Fig. Et>- The im- 
pact of the cross-correlation method on the zero-point of 
the Pp is thus significant, while the slope increases only 
slightly (in absolute value) from —0.064 to —0.08. 
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Fig. 1. Interpolated radial velocity curves based on the cross-correlation method are presented of each Ccphcid in our 
sample. Uncertainty are too small to be seen (around 0.5 km s _1 ). The horizontal lines near extrema give an indication 
about ARV C [4896] . The short horizontal lines are the 7- velocities (see Sect.[6|) corresponding to the CC method (solid 
line), the center-of-mass velocity of paper III (dotted line) and from Fernie et al. (1995, dashed line). 



We have several possible explanations for these results. 
The cross-correlation induces actually two biases: 

1. the cross-correlated radial velocities are derived us- 
ing a Gaussian fit, making the result sensitive both 
to the spectral line width (i.e. the effective tempera- 
ture and abundances) and the rotation velocity pro- 
jected on the line of sight. These two quantities, in- 
dependently, and even more the combination of both, 
are not really expected to vary linearly with the log- 
arithm of the period. It might explain why no clear 
linear relation is found between / cc and the period of 
the star. However, the mean value of the correction fac- 
tors (around 0.93 ± 0.02) have a non negligible impact 



on the zero-point of the Pp relation which decreases 
from 1.376 to 1.31 (5%). 
2. the cross-correlation method implies a mix of differ- 
ent spectral lines forming at different levels. In the Pp 
relation, the only quantity sensitive to the line depth 
is /grad (as defined in Paper II) which compares the 
amplitude of the pulsation velocity corresponding to 
the line-forming region, and the photosphere, respec- 
tively. It is thus an estimate of the velocity gradient 
within the pulsating atmosphere of the star. The Pp re- 
lation was derived in paper II for the ^896 spectral line 
which forms very close to the photosphere (D = 8%), 
while the cross-correlated radial velocity is a mix of 
thousands of spectral lines forming at different levels, 
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Table 2. The quantities used to derive the CC projection factor (p cc ) are listed. We compare also the 7-velocitics 
(Ky[cc]) derived from the CC radial velocity curves with previous results. 



Name 


V (a) 


ARV CC 


fee 


Pcc 


y 7 [GCD](b) 


Vy [N08] (c) 


V 1 [cc] 






[km s" 1 ] 






[km s _1 ] 


[km s _1 ] 


[km s _1 ] 


R TrA 


1.35±o.o3 


30.6±o.6 


0.93 ±0 .o3 


1.25±o.o7 


— 13.2±2.o 


-H.3±0.3 


— 12.2±o.6 


S Cru 


1.34±o.o3 


35.9±o.6 


0.93 ±0 .o3 


1.25±o.o7 


— 7.1±2.0 


— 3.5±o.4 


-4.7±o.6 


YSgr 


1.33±o.o3 


36.0±o.s 


0.95 ±0 .o3 


1.26±o.o5 


-2.5± 2 .o 


-l-5±0.2 


-2.4±o.6 


P Dor 


l-32±o.o2 


33.5 ±0 .2 


0.94 ±0 .o2 


l-24 ±0 .05 


7.4 ±2 .o 


9.8 ±0 .i 


8.7±o.2 


C Gem 


1-31±0.02 


27.2±o.2 


0.94 ±0 .o3 


1.23±o.o5 


6.9±2.o 


7.1± .i 


6.4±o. 2 


RZ Vel 


1.30±o.o2 


52.5±o.7 


0.91 ±0 .o2 


1-18±0.05 


24.1 ±2 .o 


24.6 ±0 .4 


24.4 ±0 .7 


i Car 


l-28±o.o2 


34.3 ±0 .2 


0.96 ±0 .o2 


1.22±o.o4 


3.6 ±2 .o 


4.4 ±0 .i 


3.6±o.2 


RS Pup 


1.27±o.o2 


46.9± .5 


0.90 ±0 .o2 


1-16±0.04 


22.1±2.o 


25.7 ±0 .2 


25.8 ±0 .5 



a The projection factor as derived in paper II. 

b The 7- velocities derived from the Galactic Cepheid Database (Fernie et al. 1995). 
c The 7-velocities or calibrated center-of-mass velocities of the stars from paper III. 



with an average depth of around D ~ 25%. The cross- 
correlation method is thus more sensitive to the veloc- 
ity gradient (because the average line depth is large), 
which may explain the increase (in absolute value) of 
the slope from —0.064 to —0.08. Moreover, in paper II 
we provided a very rough estimate of the Pp relation 
associated to the cross-correlation method considering 
only the impact of the velocity gradient (which means 
discarding the bias related to the Gaussian fit). We 
found p = [-0.075 ± 0.031] logP + [1.366 ± 0.036] (see 
Sect. 7 of Paper II). The slope we find here (-0.08) is 
consistent with this previous rough estimate of —0.075. 

These results are important to take into account when 
deriving the distance of Galactic or LMC/SMC Cepheids 
using the cross-correlation method. We emphasize that 
our Pp cc is consistent with the result by Merand et al. 
(2005), who found p = 1.27 for S Cep (P = 5.36). 

5. Wavelength dependence of the projection 
factor 

With the data at hand, we check for a possible depen- 
dence of the projection factor on the wavelength range 
used for the cross-correlation radial velocity measurement. 
For each order, we derive the cross-correlated interpolated 
radial velocity curves, and then the corresponding ampli- 
tudes ARV CC (X). Orders 59, 68 and 72 are not considered 
due to instrumental characteristics and/or unrealistic re- 
sults. For all stars, ARV CC (X) is plotted as a function of 
the wavelength, defined as the orders' average values (Fig. 
GJi). We find linear relations between these two quantities: 

ARV CC (X) = a\\ + b x , (6) 

where a\ and 6a are listed in Tab. [3] For consistency with 
the previous section the Ai?Vc C (A) quantities have been 
slightly shifted in velocity in such a way that : 

ARV CC = a A 502.2nm + b x , (7) 



where ARV CC is derived from Tab. [5] and 502. 2nm is the 
wavelength averaged over all orders. 

We also find a relation between a\ and the logarithm 
of the period of the star : 

a x = [-0.005 ± 0.001] logP - [0.002 ± 0.001] (8) 

From these results we can make two comments. First, 
the amplitude of the cross-correlated radial velocity curves 
decreases with the wavelength. From hydrodynamical 
modelling, we know that the spectral lines form over a 
larger part of the atmosphere in the infrared compared to 
optical (Sasselov et al. 1990). This effect might a key to 
understand our result: the more extended the line forming 
regions are, the lower is the amplitude of the radial ve- 
locity curves. Second, this effect is larger for long-period 
Cepheids as compared to short-period Cepheids. A reason 
might be that the mean radius, the size of the line-forming 
regions and the velocity gradient increase with the loga- 
rithm of the period. 

In order to quantify the wavelength dependency 
of the Ppcc relation, we define two correction fac- 

, /f ARV CC (A=400nm) , f 

tOrS (jA=400nm = ARV CC anCl A=700nm — 

AW A ( flyI° 0nm) )- Wc nnd the following correcting relation 
as a function of the logarithm of the period: 



/A=4oonm = [-0.01 ± 0.01] log P + [0.99 ± 0.01] , (9) 
and 



/A=700nm= [0.02 ±0.01] log P+ [1.02 ±0.01], (10) 

The reduced x 2 are respectively 0.3 and 1.3. These 
relations are shown in Fig. fU Wc find that such corrections 
are currently irrelevant given our statistical uncertainties 
on the Ppcc relation (Eq. [5]). 
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Fig. 3. (a) Wavelength dependency of the amplitude of 
the cross-correlated radial velocity curves for each star in 
our sample. The corresponding linear relation are defined 
as : ARV CC = a\\ + b\ (b) The corresponding slopes (ax) 
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Fig. 4. Corrections to apply to the Pp cc relation (Eq. O 
in the blue (A = 400nm) and in the red (A = 700nm). 

Table 3. Coefficients of the linear relations between the 
amplitude of the radial velocity curve and the wavelength. 



Name 


o\ 


bx 


R TrA 


-0.005±o.oo2 


32.93±o.80 


S Cru 


-0.005 ±0 .ooi 


38.50±o.s8 


YSgr 


-0.005±o.oo2 


38.26±o.92 


P Dor 


-0.007±o.ooi 


37.06±o.54 


C Gem 


-u.006 ± o. oi 


30.39±o.46 


RZ Vel 


-0.010±0.002 


57.76±o.88 


I Car 


-0.008±o.ooi 


38.54±o.64 


RS Pup 


— 0.011±0.002 


52.39±o.89 



f. — „r j-i ;„ J 



6. The CC 7-velocity and the K-term of Cepheids 

Interestingly, for each Cephcid in our sample, wc found 
in paper III a linear relation between the 7-vclocitics (de- 
rived using the first moment method) of the various spec- 
tral lines and their corresponding 7-asymmetries. Using 
these linear relations, we provided a physical reference to 
derive the center-of-mass 7-velocity of the stars (V^[nos]): 
it should be zero when the 7-asymmetry is zero. These 
values are consistent with a axisymmetric rotation model 
of the Galaxy. Conversely, previous measurements of the 
7-velocities found in the literature (for e.g. Fernie et al. 
1995 : the Galactic Ccpheid Database, hereafter V^[gcd]) 
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were based on the cross-correlation method, and by us- 
ing generally only few measurements over the pulsation 
cycle. These results led to an apparent "fall" of Galactic 
Cepheids towards the Sun (compared to an axisymmetric 
rotation model of the Milky Way) with a mean velocity of 
about 2 km/s. This residual velocity shift has been dubbed 
the "K-term" , and was first estimated by Joy (1939) to be 
-3.8 km/s. 




-10 10 20 30 

7-velocity (km/s) 



Fig. 5. Ky [cc] as a function of V^[gcd] (diamond) and 
V^.[no8] (triangles). The solid and dashed line are the cor- 
responding linear interpolation respectively. 

We aim to understand why such a 2km/s mistake was 
done before. An hypothesis is that the cross-correlation 
method is biased by the dynamical structure of the at- 
mosphere of Cepheids. To verify this hypothesis, we have 
the unique opportunity to compare quantitatively and in a 
consistent way V^psros], Vy[GCD] and the 7- velocities derived 
from our HARPS cross-correlated radial velocity curves 
(hereafter V^,[cc]). The comparison is done by plotting 
V^[cc] as a function of V^[gcd] and V^,[no8] respectively 
(Fig. [5]). Data are all presented in Tab.[2]and the resulting 
linear relations arc respectively : 

y 7 [CC] = [1.03 ± 0.06]V 7 [GCD] + [0.86 ± 0.78], (11) 
and 

F 7 [CC] = [1.02 ± 0.02]y 7 [N08] - [0.99 ± 0.17], (12) 

The reduce x 2 are respectively 3.0 and 3.8. 

Several conclusions must be pointed out. The slope of 
these relations are similar and close to one, which basically 
means that there is no particular trend of the 7-velocity 
with the period of the star, or at least, it remains negligi- 
ble here. As in paper III, we find a systematic difference 
of 0.86 + 0.99 ~ 1.8km s _1 between V 7 [no8] and V^jgcd], 
which is consistent with the K-term of Cepheids. However, 
the 7-velocities we derive in this paper using the cross- 
correlation method are systematically lower by (0.9 ± 0.8) 



km s _1 than the ones found in the literature (V^igcd]), and 
they are systematically larger by 1.0 ±0.2 km s _1 than the 
calibrated center-of-mass velocities (V^,[nos]). As a conse- 
quence, the cross-correlation method alone cannot explain 
alone the K-term. The CC method is sensitive in such a 
way to the dynamical structure of Cepheids' atmosphere, 
that it is responsible for ~50% of the K-term. Something 
else is requested to explain the presence of such offsets in 
previous determinations of the gamma-velocity. It could 
be related, for instance, to the quality of observations in 
the past (Joy et al. 1939) or to the different methods used 
to derive the 7-velocity (Pont et al. 1994). 

7. Conclusions 

By comparing the amplitude of our cross-correlated radial 
velocity curves with previous results based on the first 
moment method (paper II), we derived a new Pp rela- 
tion applicable to radial velocities measured by the cross- 
correlation method. This relation is crucial for the distance 
scale calibration, and in particular to derive the distances 
of LMC and SMC Cepheids (Gieren et al. 2005a; Gieren et 
al. 2009, in preparation). We find also a slight dependence 
of the Pp relation on the wavelength. Considering our cur- 
rent uncertainties this effect is negligible, but it might be- 
come significant in the near future. The next steps are 
certainly to test the impact of the signal to noise ratio, 
the spectral resolution and the metallicity on the pro- 
jection factor. The latter point will require a large sam- 
ple of Cepheids with well-measured metallicitics. These 
studies (including this work) are fully part of the inter- 
national "Araucaria Project" whose purpose is to provide 
an improved local calibration of the extragalactic distance 
scale out to distances of a few Megaparsecs (Gieren et 
al. 2005b). Moreover, the fact that the cross-correlation 
method over-estimates the amplitude of the radial veloc- 
ity curve and under-estimates the 7-velocity (compared to 
the calibrated values presented in paper III) might have 
some implications for other kinds of pulsating stars, for 
e.g. in asteroseismology. 

Moreover, we show in paper III that the K-term of 
Cepheids vanished if one considers carefully the dynam- 
ical structure of Cepheids atmosphere. From the results 
presented in this paper, we can state that the cross- 
correlation method might not be totally responsible for the 
K-term found in the previous studies (only 50% seems to 
be a consequence of the cross-correlation method) . There 
seems to be another contribution whose nature should be 
investigated. 
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Table 4. HARPS cross-correlated radial velocities (part 1) 





R TrA 




S Cru 




Y Sgr 




Dor 




C Gem 


phase 


RVcc [km/s] 


phase 


RVcc [km/s] 


phase 


RVcc [km/s] 


phase 


RVcc [km/s] 


phase 


RVcc [km/s] 


0.09 


-27.59 ± 0.65 


0.03 


-22.53 ± 0.62 


0.12 


-17.37 ±0.65 


0.02 


1.08 ±0.18 


0.04 


-4.16 ±0.33 


0.27 


-18.40 ±0.65 


0.11 


-19.98 ±0.66 


0.14 


-16.35 ± 0.44 


0.03 


0.86 ±0.17 


0.13 


-5.29 ±0.20 


0.37 


-12.05 ±0.43 


0.18 


-16.40 ± 0.42 


0.30 


-8.43 ±0.61 


0.14 


-5.94 ±0.24 


0.15 


-4.59 ±0.17 


0.39 


-11.01 ±0.56 


0.19 


-15.99 ±0.59 


0.62 


8.22 ±0.55 


0.22 


-1.35 ±0.14 


0.23 


1.10 ±0.17 


0.55 


-1.91 ±0.48 


0.34 


-7.50 ±0.49 


0.77 


17.62 ±0.46 


0.23 


-0.48 ±0.15 


0.33 


9.12 ±0.18 


0.67 


2.39 ±0.70 


0.54 


3.55 ±0.62 


0.79 


17.40 ± 0.40 


0.33 


9.28 ±0.21 


0.35 


10.54 ±0.28 


0.78 


1.36 ±0.42 


0.60 


6.51 ±0.50 


0.84 


12.60 ± 0.43 


0.40 


15.77 ±0.21 


0.37 


11.94 ±0.20 


0.79 


0.53 ±0.61 


0.76 


12.89 ±0.77 


0.95 


-15.33 ±0.39 


0.42 


17.70 ±0.16 


0.43 


16.76 ±0.30 


0.96 


-25.47 ±0.46 


0.82 


11.56 ±0.52 


0.96 


-16.71 ±0.45 


0.44 


18.76 ±0.23 


0.44 


17.31 ±0.18 


0.98 


-26.80 ±0.63 


0.91 


-10.10 ± 1.40 






0.51 


24.50 ± 0.20 


0.46 


18.69 ±0.33 






0.97 


-20.96 ±0.51 






0.52 


24.99 ± 0.19 


0.53 


21.36 ±0.27 














0.54 


26.14 ±0.28 


0.56 


21.37 ±0.17 














0.61 


27.10 ±0.27 


0.65 


15.84 ±0.39 














0.64 


25.14 ±0.25 


0.74 


5.69 ±0.51 














0.73 


10.42 ±0.18 


0.84 


1.15 ±0.17 














0.83 


1.64 ±0.14 


0.94 


0.08 ±0.17 














0.92 


1.88 ±0.16 


0.96 


-0.25 ±0.22 



Table 5. HARPS cross-correlated radial velocities (part 2) 





RZ Vel 




i Car 




RS Pup 


phase 


RV CC [km/s] 


phase 


RVcc [km/s] 


phase 


RV CC [km/s] 


0.00 


12.01 ± 1.50 


0.01 


-14.44 ±0.14 


0.02 


2.23 ±0.41 


0.05 


-1.31 ±0.67 


0.09 


-12.98 ±0.16 


0.07 


3.97 ±0.49 


0.10 


-3.62 ±0.81 


0.29 


-2.01 ±0.25 


0.12 


6.66 ±0.41 


0.36 


18.03 ±0.42 


0.31 


-0.38 ±0.15 


0.17 


9.79 ±0.48 


0.51 


28.71 ±0.32 


0.34 


1.45 ±0.18 


0.22 


13.14 ±0.36 


0.65 


44.64 ±0.51 


0.39 


4.58 ±0.14 


0.26 


16.48 ±0.37 


0.86 


43.72 ±0.77 


0.40 


4.94 ±0.13 


0.28 


17.34 ±0.37 


0.90 


37.78 ±0.60 


0.42 


6.32 ±0.19 


0.33 


20.71 ±0.45 


0.95 


29.14 ±0.55 


0.43 


6.67 ±0.18 


0.38 


24.06 ± 0.54 






0.45 


8.12 ±0.13 


0.43 


27.32 ± 0.38 






0.46 


8.45 ±0.18 


0.83 


46.81 ±0.59 






0.48 


10.08 ± 0.22 


0.90 


32.78 ± 0.65 






0.51 


11.69 ±0.18 


0.93 


15.48 ± 0.87 






0.53 


13.27 ±0.28 


0.97 


3.70 ± 0.60 






0.54 


13.17 ±0.20 










0.56 


14.68 ±0.17 










0.57 


14.51 ±0.20 










0.60 


15.67 ±0.19 










0.61 


16.69 ±0.14 










0.65 


17.59 ±0.16 










0.67 


18.10 ±0.13 










0.70 


18.48 ±0.13 










0.75 


19.28 ±0.13 










0.78 


19.61 ±0.13 










0.81 


19.20 ±0.11 










0.84 


17.57 ±0.15 










0.87 


12.12 ±0.15 










0.89 


3.45 ±0.17 










0.92 


-5.73 ±0.13 










0.95 


-10.82 ±0.12 










0.98 


-13.41 ±0.17 







